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In a recent Letter pQ, Hou et al. (HPS) pre- 
sented numerical results for the diffusion process in two- 
dimensional dusty plasma liquids with Yukawa pair in- 
teraction [2DYL], V(r) — Q 2 exp(— r/X)/r, by solving a 
Langevin equation. The mean-squared displacement 

u r (t) = (\f(t)-f(t )\ 2 )<xt 1+a (1) 

is used to distinguish normal diffusion (a = 0) from sub- 
diffusion (a < 0) and superdiffusion (a > 0). HPS ob- 
served superdiffusion and reported a complicated non- 
monotonic dependence of a on the potential stiffness 
k = a/A, where a is the mean interparticle distance. Here 
we point out that the behavior a(n) is, in fact, regular 
and systematic, whereas the observations of Ref. pQ re- 
sulted from a comparison of different system states. 

As noted in [T], a depends on k and the coupling pa- 
rameter r — (Q 2 /47T£o) x (1/afcfiT) and finding the de- 
pendence a(/c) requires to compare states with the same 
physical coupling. This can be done by fixing, for all k, 
the value r ro1 = r/r c , where T c (n) is the crystallization 
point which is well known for k < 3 [2]. For larger n, we 
obtain T c (k = 3.5) = 2340 and T c (4) = 4500. 

We have performed detailed investigations of the de- 
pendence of a on r and k [3] and observed two different 
regimes: i) for T rcl < r^ cl = 0.35, a is monotonically 
decreasing with k, at constant r . ii) for r ro1 > Tg 61 , a 
increases monotonically with k, at constant r rcl . Around 
prci _ prci^ a j g amlos t independent of k. Fig. [I] clearly 
confirms the monotonic K-dependence of a for three fixed 
values of T rcl corresponding to the parameters shown in 
Fig. 5 of [I]. 

HPS used a different coupling parameter, T c g, which 
yields an almost constant r , for n < 3. However, 
for k > 3 it corresponds to strongly varying r rcl and 
thus to different physical situations, 0], cf. top part of 
Fig. [Tl For example, their value, r e g = 100, corresponds 
to r 1 ^ = 0.76 > r o cl , for k = 3, but to T rcl = 0.24 < rf, cl , 



for k = 4. This explains the non-monotonicity of a(n) 
reported by HPS [5] . 

Thus, we report a systematic effect of screening on su- 
perdiffusion in 2DYL based on numerical simulations. An 
increase of K supports superdiffusion for r rcl < 0.6 • Tq 6 ' 
and results in an increasing diffusion exponent in this 
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FIG. f: Bottom: Exponent a vs. « for three fixed values 
of r rcl (full lines and symbols) and r e fi (dashed lines, open 
symbols, data from Ref. Q]). Top: r rcl (fi;) corresponding to 
the values r e a used in [T]. 



range of the coupling. For higher couplings, T > 
0.6 • r™ 1 , a stronger screening has the inverse effect and 
reduces the strength of anomalous diffusion. In conclu- 
sion, we have presented numerical evidence for the exis- 
tence of a monotonic dependence of anomalous diffusion 
on screening. An explanation of this behaviour is beyond 
the present Comment and will be given elsewhere. 
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